Observation of magneto-phonon resonance of Dirac fermions in graphite 
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Coherent coupling of Dirac fermion magneto-excitons with an optical phonon is observed in 
graphite as marked magnetic- field dependent splittings and anti-crossing behavior of the two cou- 
pled modes. The sharp magneto-phonon resonance occurs in regions of the graphite sample with 
properties of superior single-layer graphene having enhanced lifetimes of Dirac fermions. The greatly 
reduced carrier broadening to values below the graphene electron-phonon coupling constant explains 
the appearance of sharp resonances that reveal a fundamental interaction of Dirac fermions. 



PACS numbers: 

The quest for ultra-high quality graphene, the new cel- 
ebrated two-dimensional (2D) electron gas pQ, is driven 
by expectations of discoveries of novel physics and appli- 
cations linked to massless Dirac fermions [2j [3] . A class 
of quantum coherent effects dictated to the rather weak 
electron-phonon interaction [4j [5] is linked to the coher- 
ent coupling of the particle-hole transitions with optical 
phonons resonantly tuned by an external magnetic field 

These resonant coupling phenomena (magneto-phonon 
resonances or MPR) occur when the energy spacing be- 
tween Landau levels (LLs) is continuously tuned to cross 
the energy of an optical phonon mode. MPRs have been 
largely explored in bulk semiconductor materials [8UTU], 
in two-dimensional semiconductor systems and in quan- 
tum dots [TT|[T2]. 

Recent theoretical studies in graphene have suggested 
that MPR leads to a rich splitting and anti-crossing phe- 
nomena of the even parity E<i g long wavelength optical 
mode due to its MPR interactions with excitations across 
Dirac fermion Landau levels (magneto-excitons) [6j [7] . 

The experimental manifestation of the quantum inter- 
actions in the MPR relies on a delicate interplay between 
the lifetimes of Dirac fermion modes and the electron- 
phonon coupling. A magnetic field tunable anomaly in 
the optical phonon of graphene has been reported [321 [14] • 
In these works the resolution of the anti-crossing that is 
the signature of the coherent coupling of electron-phonon 
modes is impeded by the relatively short lifetimes of the 
fermion excitations. The manifestation of the sharp MPR 
requires damping of the Dirac fermion excitations well 
below the weak electron-phonon coupling. This quali- 
tatively new regime of interacting Dirac Fermion modes 
has not been achieved so far. 

In this letter we demonstrate the coherent mixing of 
the E2 g phonon and the Dirac fermion magneto-exciton 
mode that is the essence of the MPR. At selected values 
of the magnetic field we observe the emergence of two dis- 
tinct modes displaying the characteristic anti-crossing be- 
havior. These observations uncover unprecedented nar- 
row linewidths of Dirac fermion magneto-exciton transi- 



tions. The Landau level transition Lorentzian widths of 
3.2 meV reported here are about three times smaller than 
those of previous investigations at similar magnetic fields 
[T4HT6] and yield damping values below the graphene 
electron-phonon coupling constant. Surprisingly, the en- 
hanced Dirac fermion lifetimes required for the observa- 
tions of the characteristic MPR electron-phonon mode 
anti-crossings are found in selected regions of graphite. 

We show below that our findings can be modeled with 
a simple coupled mode Hamiltonian describing the E 2g 
mode and sharp magneto-excitons of a graphene layer. 
The sharp Landau level transitions found in this work 
and the interpretation in terms of a single-layer graphene 
model are in line with recent experimental and theoreti- 
cal studies that suggest the existence of high-quality de- 
coupled graphene flakes in bulk graphite [TSJ [T7J [18] . In 
our work the unique areas are identified by the capabil- 
ity to scan regions of the sample within the environment 
of the magneto-optics experiment. Consistent with this 
interpretation, our spatially-resolved Raman analysis re- 
veals a large non-uniformity of the graphite sample and 
areas of the sample where no MPRs are identified (see 
Supplementary). 

The main MPR is found close to 5T and involves the 
— 1^2 and —2^1 Landau level transitions of Dirac 
fermions (Landau levels in the 7r* (tt) conduction (va- 
lence) band correspond to the positive (negative) inte- 
gers) as shown in Fig. la. Representative examples of 
spectra showing the two coupled modes cj_ and cj + are 
presented in Fig.lc and d together with the magnetic- 
field evolution of the expected mode anti-crossing behav- 
ior (Fig.l b). The simple coupled mode physics underly- 
ing the data enables algebraic extraction of parameters 
of the interacting system including the Dirac fermion life- 
times. 

Raman experiments were performed at 2K in a mag- 
netic field up to 14T in a backscattering configuration 
with a tilt angle of 6 = 20° between the magnetic field 
and the c-axis direction of Kish graphite (Toshiba ce- 
ramics) which was deposited onto a Si/Si02 substrate 
(see Fig. la). A diode-pumped solid-state laser with a 
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FIG. 1: a Schematic visualizations of the phonon and 
magneto-exciton modes involved in the magneto-phonon res- 
onance (MPR). The graphene layer is illustrated in the x-y 
plane. The magnetic-field direction in the experiment is tilted 
by — 20° from the z-axis. b Energies of two coupled normal 
modes versus magnetic field close to the MPR at around B 
= 5T. c and d Raman spectra acquired at 4.55T and 5.6T, 
respectively. The energies of the two coupled modes cj_ and 
cj+ are indicated by the red and blue dots in panel b. The 
origin of the additional peak uji is discussed later in the text. 



532nm emission line was used as the excitation source 
and focussed on the sample with a spot size of about 
80/im and a power of around 30mW. The laser spot was 
moved around on the sample to select regions with opti- 
mized MPR responses. We also performed micro-Raman 
measurements (spatial resolution of 0.5/ira) at room tem- 
perature in air on a different piece of the sample which 
exhibits large spatial variations of the 2D band (see Sup- 
plementary). 

Figure 2 displays the magnetic-field dependence of the 
Raman spectra that we obtained on two different spots 
of the sample. The E 2g phonon displays drastically dif- 
ferent behavior. In Fig. 2a the spectra are very sensitive 
to the tuning of the magnetic field, while in Fig. 2b no 
changes are observed. This highlights the non-uniformity 
of the graphite and the presence of regions that behave 
like single-layer graphene. Further support for this pic- 
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FIG. 2: a and b. Magnetic field dependence of the E2 g 
phonon spectra taken at two representative regions of the 
sample. One region displays significant MPR effect while the 
other doesn't. The smooth curves are Lorentzian fits to the 
data. c. Mode energy evolution extracted from the spectra 
displayed in a. The colors of the four branches correspond to 
those used in the Lorentzian fits. 



ture comes from the spatially-dependent analysis of the 
Raman spectra that is reported in the supplementary 
material section. Behavior similar to that in Fig. 2b has 
been reported and ascribed to the usual Bernal stacked 
graphite [14] . Below we will focus on the remarkable data 
displayed in Fig. 2a. 

The magnetic-field evolution of the mode energies of 
these spectra are shown in Fig. 2c. Previous studies have 
observed a phonon energy modulation of about 40 cm -1 
using magnetic fields up to 30 Tesla [14]. Here we found 
the optical phonon energy is modulated by more than 60 
cm -1 with fields less than 6 Tesla, indicating the presence 
of very high quality Dirac fermions with long lifetimes in 
the sample. 

The extreme sensitivity of the spectral lineshape to 
magnetic field changes seen in Fig. 2a suggests the exis- 
tence of LL transitions that come into and out of reso- 
nance with the phonon energy. It is also remarkable that 
at ~6.5 Tesla the phonon becomes extremely sharp, with 
a width of 4cm _1 . To the best of our knowledge this is 
the narrowest E 2g phonon width reported in graphene- 
related materials. 

Electric-field effect studies have shown that a sharp 
E 2g phonon is linked to the absence of resonant electron- 
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hole pair transition states that the lattice vibration can 
decay into [4j[5j[T9]. This indicates that at around 6.5 
Tesla no LL transition exists that resonantly interacts 
with the phonon. It follows that observations of a very 
sharp phonon have a nontrivial implication on the LL 
structure of the underlying electronic system. Further- 
more, the significant lineshape changes around 5T sug- 
gest that a LL transition with large degeneracy comes 
into resonance with the long wavelength optical phonon. 
Comparing with current LL studies of graphitic materi- 
als [20-22 , the transition matches well with the — l—>2 
(and —2—^1) transitions of Dirac fermions. Based on 
these considerations we proceed to model the electron 
system in the MPR anticrossing displayed in Figs. 2a and 
2c as Dirac fermions with discrete LL structures charater- 
istic of those in single-layer graphene of high perfection. 

Figure 3 shows the detailed evolution of the phonon 
spectrum. We observe that the lineshape is very asym- 
metric and the asymmetry changes side when crossing 5 
Tesla. As shown, all the asymmetric spectra can be de- 
composed into two Lorentzian peaks. The red and blue 
modes which we highlight in the figure form a paired be- 
havior at about 5 Tesla. Taking into account the 20° 
angle between the magnetic field and the graphite c-axis, 
at resonance the field perpendicular to the basal plane 
is 4.7 Tesla. The resonance condition of the —1^2 LL 
transition of Dirac fermions gives (1 + =196 
meV where Ib is the magnetic length, 196 meV is the 
E2 g phonon energy. This determines the Fermi velocity 
vp to be 1.03xl0 6 m/s, in good agreement with studies 
of Dirac fermions in graphitic materials. 

The mode energy and linewidth obtained from the 
analysis are displayed in Fig. 4. The splitting of the two 
anti-crossing branches asymptotically approaches a large 
value of 70 cm -1 at the resonance field of 5 Tesla. This 
is in contrast with results in Fig. 2b and the studies re- 
ported in Ref.|14 , in which the impact of the LL tran- 
sition is seen largely as a weak perturbation that only 
weakly renormalizes the phonon energy and broadens the 
phonon Raman peak. In the results shown in Fig. 4, the 
sharpness of the LL transition leads to the strong mixing 
of the optical phonon with a magneto-exciton, and to the 
establishment of two distinct new collective modes that 
have characteristic anti-crossing behavior. 

To quantitatively analyze our data, we employ a cou- 
pled mode theory that was proposed by Goerbig et al. 
in Ref.[7 where the two coupled normal modes are de- 
scribed as 



fa± = g 2 ^ 9 

In the equation, Epn and Eme describe the phonon and 
magneto-exciton respectively, while g is the coupling pa- 
rameter. To explicitly describe the broadening of the 
modes, we let Epn = uj — 27, Eme = Q — iT where a;, f2, 
7 and T are energies and half widths of the phonon and 
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FIG. 3: Detailed Raman spectra of the main MPR observed 
for magnetic fields ranging from 4.2T to 6.4T. The theoreti- 
cal calculation of the anti-crossing associated to the magneto- 
phonon resonance (MPR) at 5T (see text) is indicated by red 
and blue dotted lines. The spectral component that corre- 
sponds to the two coupled modes is highlighted. The spectra 
are shifted vertically so that each red and blue component 
taken at a given magnetic field has its peak aligned at that 
magnetic field. 



magneto-exciton . 

At resonance uj = huj± = cj± 
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The phonon width 7 in the absence of Landau damping 
into electron- hole pairs is 4cm _1 so 7 = 2cm~ 1 . The 
analysis using data in Fig. 4 yields the energy and broad- 
ening of the two coupled modes at resonance: fwj+ = 
(1619 - il^cm- 1 and HujL = (1549 - zl^cra -1 . These 
experimental data allow us to algebraically determine 
all the parameters in the coupled mode Hamiltonian: 
uj — 1584cm -1 , T = 26cra -1 and g = 37cm _1 . 

The value of g provides a direct measurement of 
the electron-phonon coupling strength since it is di- 
rectly linked to the dimensionless electron-phonon cou- 
pling constant A [6 by g = yJ^fiujB where Hujb = 
^ h B VF . Using the experimental result (>/2 + l)huj B = 

1619+1549 cm -l ? we find A = g 36 x 1Q -3 physically, 

this coupling strength reflects the modulation rate of the 
nearest neighbor hopping integral 70 with respect to the 
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FIG. 4: MPR induced anti-crossing behavior of the coupled 
modes, a and b show the mode energy and full-width-half- 
maximum (FWHM), respectively. The blue and red lines are 
the theoretical fits of the energy of the two branches. The 
cyan line in a is the Landau level transition —1—^2 or— 2 — >• 
1. The blue and red triangles are the experimental data. The 
inset reports the magnetic-field evolution of the uji mode. 

stretching of the carbon-carbon bond length b and we 
obtain ^ = 6.66eV/A. These values are in reasonable 
agreement with experiments on monolayer graphene [23] 
and epitaxial graphene [14]. 

The Landau level transition half width of 26cm -1 
(roughly 3 meV) is an important indication that the 
Dirac fermions residing in graphite are of a very high de- 
gree of perfection. As a comparison, in epitaxial graphene 
with a reported mobility of 250,000cm 2 /(V-s) [16 , the 
Landau level width at 5 Tesla is about lOmeV. 

The narrow width T less than g, a condition not met 
before, allows for the coherent coupling between the 
phonon and the fermion excitations leading to the anti- 
crossing phenomenon seen in our data. The observation 
of the delicate MPR coherent processes in graphite opens 
further venues for fundamental optics research on Dirac 
fermions. 

We now briefly discuss the green mode labelled uoi with 
energy and width displayed in the insets to Fig. 4. The 
behavior of this mode is quite similar to the magnetic 
oscillations observed in Ref. [14] . Interestingly, we found 
that the resonance field of ojj is slightly lower (at 4.8T) 



than that of the two anti-crossing modes (at 5T). This 
might indicate disorder-induced Fermi velocity renormal- 
ization of graphene. The existence of the the ujj mode is 
a further indication of sample no n- uniformity. 

For the magnetic field range below 4 Tesla, magneto- 
phonon resonance with interband transitions of higher LL 
indexes are expected. However, we were unable to clearly 
identify such resonances, possibly due to the smaller LL 
degeneracy as well as limitations imposed by our signal 
to noise ratio detection limit. 

Starting at about 13 Tesla, the phonon splits again in- 
dicating the occurrence of another anti-crossing (Fig. 2c). 
This resonance could be due to the 0^2 LL transi- 
tion that is expected at 14.6 Tesla or the — >• 1 tran- 
sition expected at 29.1 Tesla. Because the resonance is 
incomplete, we are unable to rule out one or the other 
unambiguously. 

Supplementary materials 

Spatially-resolved Raman analysis of the 2D 
phonon peak in Kish graphite 

Kish graphite was deposited on a Si/Si02 substrate 
using the standard micromechanical cleavage method. 
Spatially-resolved Raman spectroscopy measurements 
were acquired in ambient conditions. An argon ion laser 
emitting at 488 nm, used to excite the phonons in the 
graphite sample, was focused to a spot size of about 
0.5 jam with a power of about 20 mW. The position of 
the sample was controlled by piezoelectric stages. Visual 
feedback was provided by a CCD camera, allowing for 
precise control of the position and focus of the laser on 
the sample. Raman measurements were obtained every 
1 /im, over an area of 40 /im 2 , with an acquisition time 
of one minute. 

We focus our spatially-resolved investigations on the 
2D peak. In fact, the 2D peak of graphite is extremely 
sensitive to the number of layers present in the sample 
[24} [25] . It is one of the trademark Raman signatures 
used to distinguish between single and bilayer graphene 
[24] [25] . In monolayer graphene, the 2D peak is a sin- 
gle symmetric peak whereas in bilayer graphene the 2D 
peak can be fitted with four peaks, evidence of the four 
possible double resonance processes that occur due to 
the splitting of the valence and conduction bands in the 
vicinity of the Dirac points [24-26 . Analysis of the 2D 
peak of our sample exhibits variations from the standard 
graphite peak shown in Fig. 5a. The double peak struc- 
ture shifts from including the more intense peak on the 
higher energy side to the lower energy side. The lower 
energy peak component has been shown by Lukyanchuk, 
et al., to correspond to the single layer 2D peak position 
[27] . The nonuniformity of the 2D peak is a strong in- 
dication of the decoupling of the layers leading to areas 
of the sample with a high contribution of electronically 
decoupled layers. 

The ratio between the intensity of the peak positions of 
the two components comprising the 2D peak of graphite 
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FIG. 5: a) The characteristic 2D peak of graphite, b) An- 
other 2D peak obtained on the same piece of Kish graphite 
shown in the inset, c) A map of the ratio of intensities of 
the double peak structure of the 2D Kish graphite peak was 
obtained by integrating over energy ranges centered on the 
two components of the 2D peak. The lower energy compo- 
nent was integrated from 2704 era -1 to 2718 era -1 or 9 pixels 
of the CCD camera. An equivalent number of pixels were 
selected in the high energy range, 2727 era -1 to 2742 cm -1 . 
The spatially-resolved Raman spectra were acquired with a 
step size of 1 \im. The blue regions correspond to peaks re- 
sembling panel a and the red areas or similar to panel b. 



are plotted as a function of position in Fig. 5c. The re- 
gions in which the 2D peak corresponds to the charac- 
teristic graphite 2D spectrum (Fig. 5a) are indicated by 
the blue areas in Fig. 5c. The spectra in which the lower 
energy component dominates (Fig. 5b) are red in Fig. 5c. 
Variations of the intensity of the two peaks is a clear in- 
dication of the irregular single layer contribution to the 
Raman signal. 
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